A reformulation of the Barrabes-Israel null-shell formalism 
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We consider a situation in which two metrics are joined at a null hypersurface. It often occurs 
that the union of the two metrics gives rise to a Ricci tensor that contains a term proportional 
to a Dirac 5-function supported on the hypersurface. This singularity is associated with a thin 
distribution of matter on the hypersurface, and following Barrabes and Israel, we seek to determine 
its stress-energy tensor in terms of the geometric properties of the null hypersurface. While our 
treatment here does not deviate strongly from their previous work, it offers a simplification of the 
computational operations involved in a typical application of the formalism, and it gives rise to a 
stress-energy tensor that possesses a more recognizable phenomenology. Our reformulation of the 
null-shell formalism makes systematic use of the null generators of the singular hypersurface, which 
define a preferred flow to which the flow of matter can be compared. This construction provides 
the stress-energy tensor with a simple characterization in terms of a mass density, a mass current, 
and an isotropic pressure. Our reformulation also involves a family of freely-moving observers that 
intersect the surface layer and perform measurements on it. This construction gives operational 
meaning to the stress-energy tensor by fixing the argument of the (S-function to be proper time as 
measured by these observers. 
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I. INTRODUCTION 

Many applications of the general theory of relativity 
involve two metrics being joined at a common boundary, 
a hypersurface that partitions spacetime into two distinct 
regions. It often occurs that the union of the metrics does 
not form a smooth solution to the Einstein field equation, 
even when the metric is continuous at the hypersurface. 
Indeed, a discontinuity in the metric's transverse deriva- 
tive always produces a singularity in the Riemann ten- 
sor. The singularity, however, is sufficiently mild (it is a 
Dirac (5-function with support on the hypersurface) that 
it yields itself to a sound physical interpretation. 

When the hypersurface is timelike (or spacelike), only 
the Ricci part of the Riemann tensor is singular, and 
the singularity can be associated with a material stress- 
energy tensor supported on the hypersurface. When the 
hypersurface is null, on the other hand, both the Weyl 
and Ricci parts of the Riemann tensor are singular; here 
the singularity in the Ricci tensor can just as well be as- 
sociated with a material stress-energy tensor, while the 
singularity in the Weyl tensor is interpreted as an inde- 
pendent impulsive tidal wave. The material present on 
the hypersurface, which gives rise to a singular Ricci ten- 
sor, is often referred to as a surface layer, or as a thin 
shell. The dynamics of a timelike or null surface layer, 
whose world surface coincides with the singular hyper- 
surface, can be of considerable interest and has been the 
subject of numerous studies. 

The standard formalism to describe the physical prop- 
erties of a timelike (or spacelike) surface layer was pro- 
duced by Israel more than 35 years ago ||l| . This practical 
formalism allows for efficient computations of the shell's 
surface stress-energy tensor in terms of the discontinuity 
of the extrinsic curvature across the hypersurface. The 
Israel formalism has even been implemented as the pack- 
age GRJunction |^ within the computer algebra program 



GrTensorll ||]. 

Several extensions of the Israel formalism to the case of 
a null hypersurface were proposed in the following years 
l@' H H' 0' H Hi' but a definitive treatment had to wait 
the 1991 work of Barrabes and Israel [0, who produced 
a formulation that was much better suited to practical 
applications. Their formalism was used in several inter- 
esting contexts: to model an exploding white hole [ pT| ; 
to speculate on the possible creation of expanding uni- 
verses inside black holes ||l2|; to calculate the gravita- 
tional signal that accompanies a body's sudden change 
of multipole moments |13, 0, or its sudden deceleration 
to a complete stop ||l^, 16 ; and to calculate the detection 
of an impulsive null signal by a system of neighbouring 
test particles 
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A number of additional applica- 
tions were considered by Barrabes and Hogan in their 
extension of the original Barrabes-Israel formalism [ p^ , 
which was also extended to the case of scalar-tensor the- 
ories [20|, ^ and the Einstein-Cartan theory ||2^. The 
formalism has been implemented within the computer 
package GRJunction p3| . 

Our purpose with this paper is to present a reformula- 
tion of the Barrabes-Israel formalism for null surface lay- 
ers. With this new formulation we hope to achieve both 
a simplification of the computational steps involved in a 
typical application of the formalism, and a clarification 
of some of the interpretational aspects of the formalism. 
While the work presented here does not deviate strongly 
from the original Barrabes-Israel treatment, we believe 
that it will ease and encourage further explorations of the 
physics of null surface layers; as such, our work should 
constitute a useful addition to the existing literature. In 
effect, what we set out to do is to make the Barrabes- 
Israel formalism more user-friendly, and more amenable 
to straightforward physical interpretation. 

The first key aspect of this reformulation of the null- 
shell formalism is our systematic involvement of the null 
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generators of the singular hypersurface. The generators 
are used to define a preferred flow on the hypersurface, 
to which the flow of matter on the surface layer can be 
compared. As we shall see, this construction provides the 
surface stress-energy tensor with a simple characteriza- 
tion in terms of a mass density (which represents a flow 
of matter along the null generators), a current density 
(which represents a flow of matter across generators) , and 
an isotropic pressure. The surface stress-energy tensor 
therefore comes with a clear and meaningful phenomenol- 
ogy. 

The null generators also provide us with a preferred 
family of intrinsic coordinates on the hypersurface. These 
are adapted to generators, in the sense that two of the 
coordinates serve to label the generators, while the third 
coordinate is the generator's parameter. The reformula- 
tion makes exclusive use of this family of intrinsic coordi- 
nates. While this represents a restriction of the freedom 
granted by the original treatment of Barrabes and Israel 
(who allowed for an arbitrary system of intrinsic coordi- 
nates), we shall see that this choice greatly simplifies the 
computational operations that arise in practical applica- 
tions of the formalism. A particularly important aspect 
of this simplification is that once the intrinsic coordinates 
have been constructed, all remaining operations follow a 
simple algorithm. We admit that the explicit construc- 
tion of the preferred coordinates can be a difficult task, 
and that in some cases this might motivate a retreat to 
the original formulation. In many applications, however, 
including the ones considered in this paper, the intrinsic 
coordinates are easily obtained and the remaining oper- 
ations are entirely straightforward. 

The second key aspect of our reformulation is an at- 
tempt to give operational meaning to the (5-function 
structure of the surface layer's stress-energy tensor. In 
the original treatment, the argument of the (5- function 
is left arbitrary, and as a result the stress-energy ten- 
sor acquires an ambiguity: it is defined up to a mul- 
tiplicative factor that can vary arbitrarily over the hy- 
persurface. Our reformulation does not eliminate this 
ambiguity (which is an unavoidable feature of null hy- 
persurfaces), but it characterizes it in a physical, well- 
motivated way. We introduce, in a neighbourhood of the 
hypersurface, a family of freely-moving observers that in- 
tersect the surface layer and perform measurements on it. 
This family is arbitrary, and it is this arbitrariness that 
now generates the ambiguity associated with the surface 
stress-energy tensor. Once a specific choice of observers 
has been made, however, all traces of ambiguity disap- 
pear. The argument of the (5-function is then fixed to be 
the observer's proper time, which is set to zero at the 
hypersurface. This choice confers a robust operational 
meaning to the surface stress-energy tensor. 

Our reformulation of the Barrabes-Israel null-shell for- 
malism is presented in Sec. II of the paper; there we omit 
all derivations and present the formalism as a ready-to- 
use recipe. Derivations of the main results are contained 



the parameterization of the null generators and their im- 
pact on the physics of null surface layers. Both Sees. Ill 
and [V can be omitted during a first reading of the paper; 



they contain technical material that are not required in 
practical applications of the formalism. Finally, in Sec. ^ 
we present four simple applications of the formalism, de- 
signed to illustrate the methods and highlight some as- 
pects of the physics of null surface layers. 



II. NULL-SHELL FORMALISM 

We consider a hypersurface S that partitions space- 
time into two regions A4^. In the metric is 
when expressed in coordinates indicated, the co- 

ordinate systems on each side of the hypersurface may be 
different. We assume that the hypersurface is null, and 
our convention is such that M~ is in the past of S, and 
in its future. We assume also that the hypersur- 
face is singular, in the sense that in a suitable coordinate 
system defined in a local coordinate neighbourhood that 
includes S, the Riemann tensor is distributional and con- 
tains a term proportional to a Dirac 5-function supported 
on the hypersurface. In general, both the Weyl and Ricci 
parts of the Riemann tensor are singular. The singular- 
ity in the Ricci tensor can be thought to be produced 
by a material surface layer whose spacetime history coin- 
cides with S. We wish to give a characterization of this 
surface layer in terms of a surface stress-energy tensor. 
We shall here have nothing to say about the singularity 
in the Weyl tensor and the interesting physics associated 
with it; for a discussion we refer the reader to the relevant 
literature (for example, see [|, |, ||l0[ |l|). 



Intrinsic coordinates 

Our formulation of the null-shell formalism is based on 
a preferred family of coordinate systems that are used to 
represent events on the hypersurface; these coordinates 
are adapted to the null generators of S. Since the genera- 
tors define a preferred flow on the hypersurface, and since 
the preferred coordinates give the simplest description of 
the generators, this type of intrinsic coordinate system 
will be seen to produce an especially simple description 
of the surface stress-energy tensor. 

The preferred intrinsic coordinates are 



y 



(2.1) 



in Sec. III. In Sec. IV we examine some issues regarding 



where the lower-case latin index a runs over the values 
(1,2,3), while the upper-case index A runs over (2,3); 
thus = X, — 9"^ , and = 9^. We assume that the 
coordinates y"" are the same on both sides of S. We take 
A to be an arbitrary parameter on the null generators 
of the hypersurface, and we use 6^ to label the genera- 
tors. Thus, to each generator we assign two coordinates 
9^, and to each event on a given generator we assign a 
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third coordinate A. Changing A while keeping 6"^ con- 
stant produces a displacement along a single generator; 
changing 9^ while keeping A constant produces a dis- 
placement across generators. It is possible to choose A to 
be an affine parameter on one side of the hypersurface; 
as we shall see in Sec. IV, however, it is in general not 
possible to make A an affine para mete r on both sides of 
S. The coordinate system of Eq. (2.1) is not unique, as 
it leaves the freedom of performing an independent repa- 
rameterization of each generator; we shall explore this 
freedom in Sec. 



IV 



Tangent vectors 

As seen from , E is described by a set of parametric 
relations a:^(j/''), and using these we can introduce tan- 
gent vectors e^^ — dx^./dy"' on each side of the hyper- 
surface. These are naturally segregated into a null vector 
that is tangent to the generators, and two spacelike 
vectors that point in the directions transverse to the 
generators. Exphcitly, 



(2.2) 



(Here and below, in order to keep the notation simple, we 
refrain from using the "±" label in displayed equations; 
this should not create any confusion.) By construction, 
these vectors satisfy 



kfy^k ' ■ /Cq6^. 

The remaining inner products 



aAB{\e^) = 



a 

9a0 e^e^ 



(2.3) 



(2.4) 



do not vanish, and we assume that they are the same on 
both sides of E: 



[<Jab] - 0. 



(2.5) 



We have introduced here the standard notation [A] = 
A{M^)\j: — A{M~)\i: for any scalar quantity A. In the 
first term, A is computed in A4~^ and evaluated on E; 
in the second term, A is computed in instead. A 
vanishing [A] means that A is continuous at the hyper- 
surface; a nonzero [A] means that A is discontinuous, and 
[A] gives the jump of A across E. 



Intrinsic metric 

It is easy to see that aAB acts as a metric on E: A 
displacement on the hypersurface is given by dxj. = 

we find 



and by virtue of Eqs. (2.3) and (2.4) 



ds'i. 



aAB dO'^dO^ 



(2.6) 



on either side of S. The condition ( |2.5| ) therefore en- 
sures that the hypersurface possesses a well-defined in- 
trinsic geometry, a central requirement of the null-shell 
formalism. A violation of Eq. (2.5) would indicate that 
the singularity of the hypersurface is of a much stronger 
type than can be handled by distributional techniques, 
and the physical interpretation of such a singular hyper- 
surface would be much more delicate. 

It should be noted that the intrinsic metric of the null 
hypersurface is two-dimensional. This is as it should 
be, as a displacement along a generator necessarily gives 
ds^ — 0. This property, that the hypersurface metric is 
not merely degenerate but explicitly two-dimensional, is 
a definite advantage. For example, there is no difficulty 
in defining a metric inverse, which we denote o"^^, and 
which is the same on both sides of the hypersurface. The 



inverse two-metric satisfies a^'^ 



"c- 



Transverse vector 

The tangent vectors fc^ and e^^ provide only a par- 
tial basis for the decomposition of vector fields on either 
side of the hypersurface. We now complete the basis by 
adding a transverse null vector that satisfies 



NaN" = 0, 



-1, 



0. (2.7) 



These four equations give a unique specification of the 
transverse vector. The basis gives us the completeness 
relations 



,AB, 



A'=-B 



(2.8) 



for the inverse metric on either side of E (in the coor- 
dinates x%)\ as was noted above, a^^ is the inverse of 
uab- The completeness relations can be used to decom- 
pose a vector A'^ in the specified basis: A" = g'^'^Ap = 

{-N^A)3)k" + {-kPAp)^ + (o-^^e§A^)e^, where we 
have omitted the ± label for increased clarity. 



Congruence of timelike geodesies 

The first key aspect of our reformulation of the null- 
shell formalism was introduced previously: it is our 
choice of preferred intrinsic coordinates which, as we shall 
see, provides the surface stress-energy tensor with a sim- 
ple description and a clear physical interpretation. The 
second key aspect of our reformulation is the introduction 
of an arbitrary congruence of timelike geodesies that in- 
tersect the null hypersurface. As we shall see, the congru- 
ence will give operational meaning to the distributional 
character of the material stress-energy tensor. 

Each member of the congruence is thought of as the 
world line of a freely-moving observer that intersects the 
surface layer and performs measurements on it. The con- 
gruence then corresponds to a whole family of observers. 
We treat all such families equally and do not attempt to 
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introduce a preferred set of observers; the congruence is 
therefore completely arbitrary. 

The timelike geodesies are parameterized by proper 
time T, which is adjusted so that t — when a geodesic 
intersects E; thus, r < in Ai~ and r > in 
The vector field tangent to the congruence is u^. and 
is normalized; a displacement along a single geodesic is 
given by dxj. = dr. To ensure that the congruence is 
smooth at the hypersurface, we demand that be "the 
same" on both sides of E. This means that UaC", the 
tangential projections of the vector field, must be equal 
when evaluated on either side of the hypersurface: 



(2.9) 



If, for example, u" is specified in A4 , then the three 



conditions of Eq. (2.9) are sufficient (together with the 
geodesic equation) to determine the three independent 
components of in M+. We note that -UaN°', the 
transverse projection of the vector field, is allowed to be 
discontinuous across E. 

The proper-time parameter on the timelike geodesies 
can be viewed scalar field t{x'^) defined in a 

normal convex neighbourhood of E: Select a point a;^ 
off the hypersurface and locate the unique member of 
the congruence that links this point to E; the value of the 
scalar field at Xj. is equal to the proper-time parameter of 
this geodesic at that point. The hypersurface E can then 
be described by the statement T{x'j.) — 0, and its normal 
vector will be proportional to the gradient of t{x'^) 
evaluated on E. It is easy to check that the appropriate 
expression is 



dr 
dx^ 



(2.10) 



We recall that the factor —kf^u^^ in Eq. (2.1C) is contin- 
uous at E. 



Transverse curvature 

We have seen that the two-dimensional induced metric 
(TAB is necessarily the same on both sides of the hyper- 
surface. In the presence of a surface layer, however, we 
can expect that the metric's transverse derivative will be 
discontinuous across E. In Sec. [II we show that informa- 



tion about the metric's transverse derivative is properly 
encoded in the transverse curvature C^i defined by 



C'ab — ~Na e^.aejj — Cba- 



(2.11) 



In the presence of a surface layer, this object (which 
transforms as a scalar under transformations of the co- 
ordinates Xj.) is discontinuous across E: [Cab] 7^ 0. 
The symmetry of C^^ follows from the identity e^.^e^ = 
e".^ef , which states that each basis vector e" is Lie trans- 
ported along any other basis vector. Note that Cab is the 
same object that was denoted ICab in the original work 
by Barrabes and Israel pOl]. 



Surface stress-energy tensor 

The jump in the transverse curvature is directly related 
to the stress-energy tensor of the surface layer. To display 
this relationship we introduce the surface quantities 



AB 



AB 



AB 



XB 



(2.12) 
(2.13) 
(2.14) 



We prove in Sec. |I| that the layer's surface stress-energy 
tensor is given on cither side of the hypersurface by 



S'"'3 = fik^k^^ + j^(fc"e^ + e'Xk'^) + 



pa^^e1e%. (2.15) 



In this expression, the first term represents a flow of mat- 
ter along the null generators of the hypersurface, and the 
surface quantity therefore represents a mass density. 
The second term represents a flow of matter in the direc- 
tions transverse to the generators, and the surface quanti- 
ties therefore represent a current density. Finally, the 
surface quantity p clearly represents an isotropic pres- 
sure. The surface layer can therefore be characterized 
by four quantities: a density /i, a current j"*, and an 
isotropic pressure p. This provides the surface stress- 
energy tensor with a transparent physical interpretation. 
As was noted before, this clarity of interpretation comes 
from employing the null generators to define a preferred 
flow on the hypersurface, to which the flow of matter can 
be compared. 



We show in Sec. HI that the complete stress-energy 
tensor of the surface layer is given by 



(2.16) 



In this expression, the factor (— /c^u^)^^ is continuous at 
E, and (as was discussed previously) the vector field is 
tangent to an arbitrary congruence of timelike geodesies 
parameterized by r. The singular character of the stress- 
energy tensor is revealed by the Dirac 5-function, and 
the congruence of timelike geodesies was introduced for 
the specific purpose of making the argument of the S- 
function meaningful. Indeed, the fact that it is r — 
proper time as measured by a family of freely-moving 
observers intersecting the hypersurface — and nothing 
else confers immediate operational meaning to the pref- 
actor {—kfiU^)^^S'^. This choice of argument therefore 
clarifies the physical meaning of the surface stress-energy 
tensor. 

It is important to note that the presence of the fac- 



tor 



-k^,u^')~ 



in front of S'^ in Eq. (2.16) implies that 



/i, j , and p are not truly the surface quantities that 
would be measured by the freely- moving observers. The 
physically-measured surface quantities are instead given 

by 



(2.17) 
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1 -A 



Ppm = {-kf^u'^) V- 



(2.18) 
(2.19) 



It is remarkable that the ambiguity associated with the 
arbitrary choice of congruence is hmited to a single mul- 
tiplicative factor. The "bare" quantities /i, j^, and p are 
independent of this choice, and it is often more conve- 
nient to work in terms of those. 



User manual 

Our reformulation of the Barrabes- Israel null-shell for- 
malism is now complete. The formalism can be used 
as follows. First, locate the null hypersurface S in 
and construct the intrinsic coordinates y"" = (A, 9^); this 
step is the only nonalgorithmic operation involved in the 
recipe. Second, construct the tangent vectors A:^, e^^ 
and complete the basis by computing the transverse null 
vector . Third, compute the intrinsic metric ctab and 
verify that it is the same on both sides of the hypersur- 
face. Fourth, compute the transverse curvature and 
its jump across S. Fifth, and finally, compute the surface 
quantities /i, j^, and p. 



III. DERIVATIONS 

In this section we provide derivatio ns of t he ma in re- 
sults presented in Sec. II, namely Eqs. ( 2.12 )-( ^J^). Fol- 
lowing the Appendix of the Barrabes-Israel paper [[10[ , we 
introduce a suitable coordinate system in a coordinate 
neighbourhood that includes the null hypersurface S and 
extends into both regions A4^. These coordinates, which 
are distinct from the systems a;^ encountered before, will 
be used for the derivations, but as we saw in Sec. II, the 
final formulation of the null-shell formalism will not refer 
to them. 



the metric might be discontinuous at the hype rsurf ace. 
To compute the Riemann tensor we use Eq. (2.10) to 
relate the gradient of r to the null vector k" . The result 
is 



e(r) + i?_V0(-^) + ^s 



(3.2) 



where R^^p^s denotes the Riemann tensor in con- 



structed from g^a, and 



R 



/3S\ kj 



(3.3) 

is the singular part of the Riemann tensor. Here, [F"^^] 
denotes the jump in the affine connection across S. 
We must now characterize the discontinuous behaviour 

of 

5a/3,7- The condition [gafs] — guarantees that 
the tangential derivatives of the metric are continuous: 
[gaf3,j]k'' = = [^0/3,7] Cq. The only possible discontinu- 
ity is therefore in gafs^-yN'^ , the transverse derivative of 
the metric. In view of Eq. (2.7) we conclude that there 
exists a tensor field 7^/3 such that 



-la/3k^ 



(3.4) 



This tensor is given explicitly by "/af3 = \ga(i,-^^'^ and it 
follows that 



[^7] = -2 + l%kf, - 7^^fc«). (3.5) 

Substituting this into Eq. ( ^.3[ ) gives 

^s"/375 = ^(-fcM"^)"^(7"5fc/3fc7 -7/35^"fc7 

-7"fc/3fc5 + 7/37fc"^5)'5M: (3-6) 



and we see that fc" and ^^13 give a complete characteri- 
zation of the singular part of the Riemann tensor. 



Surface stress-energy tensor 



Singular part of the Riemann tensor 

In the domain of our suitable coordinate system we 
write the metric as a distribution-valued tensor: 



9afi = 5a/30W + 5a/30(-^)' 



(3.1) 



where 6(t) is the Heaviside step function and g^p{x^^ 
is the metric in A^='=. To ensure that the metric is well 
defined at t = we assume that in the coordinates . 



[Sa/?] = 0; Eq. (2.5) is compatible with this requirement. 
We also impose [fc^] = [e^] = [TV"] = [u"] = 0, and these 
continuity statements provide a precise meaning to the 
phrase "suitable coordinate system." 



The metric of Eq. (3.1) generally gives rise to a singular 



Riemann tensor, even when continuity of g^^ at t = is 
enforced; the reason is that the transverse derivative of 



From Eq. (3.6) we form the singular part of the Ein- 
stein tensor, and the Einstein field equations give us the 
singular part of the stress-energy tensor: 



(3.7) 



where 



S' 



a/3 



IGtt 



{k°'-i^^k^' + kP-^°'^k^' -Y^,k°'k^ -It.uk^'k" g'^P) 

(3.8) 

is the surface stress-energy tensor of the null shell, up to 
a factor {—k^u^) that depends on the choice of observers 
making measurements on the shell. 

The expression for the surface stress-energy tensor 
can be simplified if we decompose it into the basis 
(/c", A^", ej^). For this purpose we introduce the projec- 
tions 



7A = 7a/3 elkl^, 



lAB = 70/36^63: 



(3.9) 
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and we use the completeness relations of Eq. (^^) to find 
that the vector 'y°'^k^ admits the decomposition 



2V r p 



(3.10) 



Substituting this into Eq. (3.8) and involving once more 



the completeness relations, we arrive at our final expres- 
sion for the surface stress-energy tensor: 



where 



{k'^e'\ + e\k^) + 



pa 



1 

16^ 



is the surface density, 

1 

16^ 

is the surface current, and 



■A _ 



J 



'o^^^ab) 



■ AB \ 
cr 7B 1 



V 



(3.11) 



(3.12) 



(3.13) 



(3.14) 



is the isotropic surface pressure. 

The surface stress-energy tensor of Eq. (3.11) is ex- 



pressed in the suitable coordinates . This is, however, 
a tensorial equation involving vectors (A;" and e^) and 
scalars (/i, j'^, and p). This equation can therefore be 
expressed in any coordinate system. In particular, when 
viewed from A^^, the surface stress-energy tensor can 
be expres sed i n the original coordinates cc^. This means 
that Eg. (3.7) is formally ide ntica l to Eq. ( ^.16 ), while 
Eq. ( |3.1lD is identical to Eq. (|2.15| ). 



Intrinsic formulation 



The surface quantities of Eqs. ( [3.12 )-(3.14) are con- 
structed from the tensor 7^^ which is defined only in the 
suitable coordinates x'^ . As it might be grossly impracti- 
cal to construct such a coordinate system, it is important 
to find an independent way of obtaining these quantities. 
It is at this stage that we introduce the transverse curva- 
ture Cah defined by Eq. ( |2.11| ). Recalling from Eq. ( pl^ ) 
that NaC" are constants (cit 



ler zero or minus one), we 



express the transverse curvature as 



/3 



Cab = A^Q;/3e"eJ; 



(3.15) 



and we seek to relate its jump to the surface quantities. 
In the next paragraph we will use the suitable coordinates 
x^ to establish the general validity of the scalar relations 



(2.12)-(2.14). This will bring our derivations to a close. 

In the suitable coordinates, the jump in the transverse 
curvature is given by 



(3.16) 



where we have used Eq. (3.5) and the fact that fc" is 
orthogonal to e". We therefore have [C\\\ = 



2-7a/3 

[Ca\] = ^lape'Xk'^ = \lA, and [Cab] = ^lape'^e^ = 
hlAB : wh er e we have involved Eq. ( |3.9[ ). Finally, using 



Eqs. ( 3.12| )-(3.14) we find that t he su rf ace q uantities can 
indeed be expressed as in Eqs. ( 2.12 )^( 2.14 ). 

We have established that the shell's surface quanti- 
ties can all be computed directly in terms of the induced 
metric aAB and the jump of the transverse curvature Cab 
across S. This confirms that the suitable coordinates x^ 
are not needed in practical applications of the null-shell 
formalism. 



IV. PARAMETERIZATION OF THE NULL 
GENERATORS 

In this section we discuss two technical issues related to 
the parameterization of the null generators of the singular 
hypersurface. In the first subsection we justify a state- 
ment made earlier, that in general A cannot be an afhne 
parameter on both sides of the hypersurface; our presen- 
tation here is based directly on Sec. Ill of the Barrabes- 
Israel paper [ p^ , but for completeness we adapt their 
discussion to the notation of this paper. In the second 
subsection we examine the transformation properties of 
the surface quantities under a reparameterization of the 
null generators. 



Afflne parameterization of the null generators 

Whether or not A is an afline parameter can be decided 
by computing k±, the "acceleration" of the null vector 
/c^. This is defined on either side of the hypersurface by 



(4.1) 



and A will be an affine parameter on the M.^ side of S 
if K± = 0. According to Eq. dJ), k = -Nak^'.^k'^ = 

—Ng e^.ge ^ — C\x, where we have also used Eqs. (2.2) 



and (2.11). Equation (2.14) then relates the discontinuity 
in the acceleration to the surface pressure: 



in p. 



(4.2) 



We conclude that A can be an affine parameter on both 
sides of S only when the surface layer has a vanishing 
surface pressure. When p 7^ 0, on the other hand, A can 
be chosen to be an affine parameter on one side of the 
hypersurface, but it will not be an affine parameter on 
the other side. 

Additional insight into this matter can be gained from 
Raychaudhuri's equation. 



d9 



1 



(4.3) 



which describes the evolution of the congruence of null 
generators [Bil l2^; the equation holds on either side of 



7 



the hypersurface. Here, 9 and ctq^ are the expansion and 
shear of the congruence, respectively, and we have set 
the vorticity to zero because the congruence is hypersur- 
face orthogonal. Because it depends only on the intrinsic 
geometry of the hypersurface, the left-hand side of Ray- 
chaudhuri's equation is necessarily continuous across the 
shell. Continuity of the right-hand side therefore implies 



]6i = 87r[T„/3fc"/c^]. 



(4.4) 



This relation shows that [k] ^ (and therefore p ^ 0) 
whenever the component Tapk'^k^ of the stress-energy 
tensor is discontinuous across the shell. We conclude 
that A cannot be an afRne parameter on both sides of S 
when [Tapk^k^] ^ 0. (Notice that this conclusion breaks 
down when 6* = 0, that is, when the shell is stationary.) 

Recalling that the expansion 6 is equal to the fractional 
rate of change of the congruence's cross-sec tion al area, 
we find that with the help of Eq. (U), Eq. ( |4^ ) can be 
expressed as —p9 = [Ta^/c^fc'^], or 



where fc" = ((9a;«/9A)0^ and = {dx'^/de^)x, or by 

dx°' ^k°'dX + e'Xde^, (4.10) 

where fc" = {dx°'/dX)g^ and = {dx" /d9^))^; these 
relations hold on either side o f S , in the relevant coordi- 
nate system x^. Using Eq. ([f.7|), it is easy to see that 
the tangent vectors transform as 



-Pk° 



eA - CAB 



(4.11) 



under the reparameterization of Eq. ( |4.6| ). It may be 
checked that the ne w ba sis vectors satisfy the orthogonal- 
ity relations of Eq. (2.3), and that the induced metric is 



invariant under this transformation: aAB = gap^A^s 
.gct/3e^e^ = ctab- To preserve the relations (2.7) we let 
the new transverse null vector be 



1 



c'^CAC-'^ k" 



(4.12) 



p—dS ^[T^0k''k^]dS, 



(4.5) 



where dS is an element of area transverse to the genera- 
tors. This equation has a straightforward interpretation: 
The left-hand side represents the work done on the shell 
as it expands or contracts, while the right-hand side is 
the energy absorbed by the shell from its surroundings; 
Eq. (4.5) therefore states that all of the absorbed energy 
goes into work. 



Reparameterization of the null generators 

We now investigate how a change of parameterization 
affects the surface density /z, surface current j^, and sur- 
face pressure p of the null shell. Because each generator 
can be reparameterized independently of any other gen- 
erator, we must consider transformations of the form 



A ^ A(A,6i^). 



(4.6) 



The question before us is: How do the surface quantities 
change under such a transformation? 

To answer this we need to work out how the transfor- 
mation of Eq. ( [4.6| ) affects the vectors fc", e" , an d N°'. 
We first note that the differential form of Eq. (|4.6|) is 



d\ ^el^ d\ + CA d9^ 



where 



dX 



CA 



{d9^\' 



(4.7) 



(4.8) 



both and ca depend on = (A, 9^), but because they 
depend on the intrinsic coordinates only, we have that 
[e^] = = [ca]- a displacement within the hypersurface 
can then be described either by 



dx" =k°'d\ 



e'Xd9^, 



(4.9) 



where c = a cb - This ensures that the completeness 
relations of Eq. ( |2.8| ) take the same form in the new basis. 

It is a straightforward (but slightly tedious) task to 
compute how the transverse curvature Cab changes un- 
der a reparameterization of the generators, and to then 
compute how the surface quantities transform. We find 



-A 



e/fi + 2cAj^ + c^CAe-^p, 



e-f'p. 



(4.13) 
(4.14) 
(4.15) 



These transformations, together with Eq. (4.11), imply 
that the surface stress-energy tensor becomes S°'^ — 
e-t^S"'^. We also have (-A:^u'^)-i = e''(-fc^M^), and 
these results imply that the combination {—ki^iu'^)"^ S"^ 
is invariant under the reparameterization. This, finally. 



establishes the invariance of T, 



Q/3 



As a final rema rk, we note that under the reparame- 
terization of Eq. (|4.6|), t he ph ysically-measured surface 
quantities of Eqs. (2.17)-(2.19) transform as 



.2/3 



A 

C CAPpm 1 



(4.16) 
(4.17) 
(4.18) 



we see in particular that the physically-measured surface 
pressure is an invariant. 



V. APPLICATIONS 

In this section we consider four simple applications of 
the null-shell formalism. Those are designed to illustrate 
the basic methods involved and highlight some aspects 
of the physics of null surface layers. 
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Imploding spherical shell 

For our first application of the null-shell formalism, 
we look at the gravitational collapse of a thin spherical 
shell. We imagine that the collapse proceeds at the speed 
of light, so that the world surface of the shell coincides 
with a null hypersurface S. We take spacetime to be flat 
inside the shell (in M.^), and write the metric there as 



(5.1) 



in terms of spatial coordinates (r, 9, (j)) and a time coor- 
dinate t-. The metric outside the shell (in A^+) is the 
Schwarzschild solution. 



dsi 



-fdtl + f-'^dr^+r^ dQ^ 



(5.2) 



expressed in the same spatial coordinates but in terms 
of a distinct time we have / = 1 — 2M/r and Af 
designates the gravitational mass of the imploding shell. 

As seen from Ai^ , the null hypersurface is described 
by the equation t_ + r = v_ — constant, which implies 
that the induced metric on E is given by ds^ — dil^. 
As see from A^+, on the other hand, the hypersurface is 
described by t+ -I- r*(r) = v+ = constant, where r*(r) = 
J /-I dr = r + 2M ln(r/2M - 1), and this gives rise to 
the same induced metric. From these considerations we 
conclude that it was permissible to express the metrics 
of A^* in terms of the same spatial coordinates (r, 6, cj)), 
but that cannot be equal to t_. The shell's intrinsic 
metric is 



(TAB d9^d9^ = (de^ + sin^ edc/)^), 



(5.3) 



where we have set 6^ = [9, cf)) and identified —r with the 
parameter A on the null generators of the hypersurface. 
We shall see that here, A is an afhne parameter on both 
sides of E. 

As seen from A^~, the parametric equations = 
a;" (A, 0^) that describe the hypersurface have the explicit 
form i_ = w_ + A, r = —A, 9 — 6, and (p ~ cf). These 
give us the tangent vectors k"da = dt — dr, e^da = do, 
and e'^da = d^, and the basis is completed by the trans- 
verse vector Na dx" — ~^{dt ~ dr) . From all this and 
Eq. ( 2.11 ) we find that the nonvanishing components of 
the transverse curvature are 



(5.4) 



The fact that C^^ — confirms that A = — r is an affine 
parameter on the side of E. 

As seen from + , the parametric equations arc — 
r*(— A), r = —A, 9 = 6, and (p = (p. The basis vectors 
are k'^da = f^^dt - dr, Cgda = dg, e'^da = d^, and 
N^dx" = ~^{f dt ~ dr). The nonvanishing components 
of the transverse curvature are now 



2r 



C^AB- 



(5.5) 



The fact that C^-^ = confirms that A = — r is an affine 
parameter on the A4~^ side of E; A is therefore an affine 
parameter on both sides. 

The angular components of the transverse curva- 
ture are discontinuous across the shell : [Ca b] = 
-{M/r^)aAB- According to Eqs. (|2.12| )- (|2.14|) , this 



means that the shell has a vanishing surface current 
and a vanishing surface pressure p, but that its surface 
density is 



M 

47]- 7-2 



(5.6) 



We have therefore obtained the very sensible result that 
the surface density of an imploding null shell is equal 
to its gravitational mass divided by its (ever decreasing) 
surface area. Notice that /ipm = A* for observers at rest in 
Ai~ . Because of the gravitational action of the null shell, 
however, these observers do not remain at rest after cross- 
ing over to the A4~^ side: A simple calculation, based on 
Eq. (p!s|), reveals that an observer at rest before crossing 
the shell will move according to dr/dr = —(7^ — /)^^^ af- 
ter crossing the shell; the constant 7 varies from observer 
to observer, and is related by 7 = 1 — M/r-^ to the radius 
rj] at which a given observer crosses the hypersurface. 



Spherical null shells 

In our second application of the null-shell formalism, 
we extend the results of the preceding subsection to a 
situation in which the surface layer is immersed in a gen- 
eral spherically-symmetric spacetime. Our results here 
are identical to those first presented by Barrabes and 
Israel in Sec. IV of their paper |jl^ . We adapt their dis- 
cussion to the notation of this paper for completeness, as 
spherically-symmetric situations are of extreme practical 
importance. We assume that the surface layer is either 
expanding or contracting, and we exclude stationary hy- 
persurfaces from our considerations. For a discussion of 
the interesting special case of horizon-straggling shells, 
we refer the reader to the original Barrabes-Israel work. 

The metric in Ai^ is expressed in the general form 



3'^'dw{fe'''dw + 2Cdr)+r^dn^, (5.7) 



where ip± and f± are two arbitrary functions of the co- 
ordinates w± and r; we also have the parameter C = ±1 
whose meaning will be explained shortly. (Notice that 
the same r is used as a coordinate in both regions.) It 
is convenient to introduce the mass functions m± {w± , r) 
defined by the relation 



/ = 1 - 2m/r 



(5.8) 



The hypersurfaces w± = constant are null. They are 
expanding (r increases along the generators) if C = +1, 
and w± (which normally would be denoted u±) is then 
properly interpreted as a retarded-time coordinate. On 
the other hand, the null hypersurfaces are contracting 
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(r decreases along the generators) if C = — and w± 
(which normaUy would be denoted v±) is then properly 
interpreted as an advanced-time coordinate. 

The boundary between A4~ and A^"*" is taken to be 
a pa rticu lar hypersurface w± = constant. According to 
Eq. (5.7), its intrinsic metric is given by 



aAB de^de^ = {de^ + sin^ 6 d(l)^ 



(5.9) 



which is the same on both sides of S. We have identified 
9^ with the two angles {9,(f>), and we take X = (r to 
be the parameter on the null generators. Notice that A 
always increases along the generators. 

The parametric equations that describe the hypersur- 
face are w± = constant, r = (X, 9 — 9, and (f) — (f). The 
tangent vectors are then 

fc"5„=C5., e^do^^de, e%d^ = (5.10) 

on either side of the hypersurface. The basis is completed 
by the transverse null vector 



No, dx" - i {fe^' dw + 2C dr) . 



(5.11) 



We may now use Eq. (2.11) to compute the transverse 
curvature. The nonvanishing components are 



dr ' 



CAB = ~C^aAB- (5.12) 

Zr 



We recall from Sec. EV that k± is the "acceleration" of 



the null vector A:^; from the first equation we therefore 
learn that A is an affine parameter on the side of E 
when d'ip±/dr = 0. 

According to Eqs. (|t|)-(|t|) and (|5.12|) , the only 
nonvanishing surface quantities are /x and p. Using 
Eq. (5^) we find that the mass density is given by 



= (-0 



47^7-2 ' 



(5.13) 



this is the difference in the mass function across the shell, 
[m](r), divided by the shell's surface area, 47rr^. Assum- 
ing that /J, is positive, we see that [m] < if the hyper- 
surface is expanding — the shell removes energy from 
the centre — whereas [m] > if the hypersurface is con- 
tracting — the shell then brings energy to the centre. 
The shell's surface pressure is given by 



dr 



(5.14) 



This result can be related to the first law of thermody- 
namics if we note that [dtp/dr] = ^■Kr\Taj3k°'k^] accord- 
ing to the Einstein field equations, and that the expan- 
sion of the null generators is given by — k^.^ — n = 
Equation (5.13) is then seen to be equivalent to the state- 
ment p9 — —[Tapk^k^] previously displayed in Eq. (4.5) 



Accreting black hole 

In our third application of the null-shell formalism, we 
consider a nonrotating black hole of mass (M — m) which 
suddenly acquires additional material of mass m and an- 
gular momentum J = aM <C M'^. We describe the ac- 
creting material by a singular stress-energy tensor sup- 
ported on a null hypersurface E. 

The spacetime in the future of S — in — is that 
of a slowly rotating black hole of mass M and (small) 
angular momentum aM. We write the metric in A4^ as 



dsl 



-fdt^ + f~^dr^ 



■dn^ 



AMa 



' dtd<t), 



(5.15) 

where / = 1 — 2M/r\ this is the slow-rotation limit of 
the Kerr metric, and throughout this subsection we will 
work consistently to first order in the small parameter a. 

As seen from A^+, the null hypersurface is described by 
V = t + r* = 0, where r* = J dr = r + 2M\n{r/2M- 
1). In the slow-rotation limit, every surface v = constant 
is null, and it follows that the vector /c" — g°'^{—dpv) is 
normal to E and tangent to its null generators. We have 



k^d^ 



dt 



Or + ^rrd^, 



(5.16) 



and from this expression we deduce four important prop- 
erties of the generators. First, the generators are affinely 
parameterized by A = —r. Second, as measured by in- 
ertial observers at infinity, the generators move with an 
(ever increasing) angular velocity 



dt 



2Ma 



generators 



(5.17) 



Third, 9 is constant on each generator. And fourth, in- 
tegration of d(p/{—dr) — 2Ma/{r^ f) reveals that 



V' = 0+-(i 

r 



(5.18) 



also is constant on the generators. 

We shall use = (A = —r,9,ip) as coordinates on E, 
as these are well adapted to the generators. Remember- 
ing that dt = -dr/ f and d4> = dip — {2Ma/r^f) dr on E, 
we find that the induced metric is 



aAB d9'^d9^ = {d9^ + sin^ 9di!^), 



(5.19) 



and we see that the hypersurface is intrinsically spherical. 

The parametric description of E, as seen from A^ + , is 
a;"(— r, 0, ?/)), and from this we form the tangent vectors 
e" — k", eg — Sg, and ej^ — ^4,- The basis is completed 
by the tra nsver se null vector N^dx" — ^{~fdt + dr). 
From Eq. (2.11) we obtain 



C 



Mia 



sin^ 9, 



C 



AB 



2r 



(^AB 



(5.20) 



for the nonvanishing components of the transverse cur- 
vature. 



10 



The spacetime in the past of S — in — is that 
of a nonrotating black hole of mass (M — m). Here we 
write the metric as 



^d^P^), (5.21) 



in terms of a distinct time coordinate t and the angles 9 
and ■0; we also have F = 1 — 2{M — m)/r. This choice of 
angular coordinates reflects the fact that zero-angular- 
momentum observers in corotate with the shell's 
null generators; this is a manifestation of the dragging 
of inertial frames. As we shall see presently, this choice 
of coordinates is dictated by continuity of the induced 
metric at S. 

The mathematical description of the hypersurface, as 
seen from Ai~ , is identical to its external description 
provided that we make the substitutions t ^ t, cj) ^ iIj, 
M ^ M — m, and a — > 0. According to this, the induced 
metric on E is still given by ds^ = r^(c?6'^-f sin^ 9 dijP'), as 
required. The basis vectors are now k^da — F~^^dt — dr, 
e^da = de, e'^da = d^, and Nadx" = ^{—Fdi + dr). 
This gives us 



(5.22) 



for the nonvanishing components of the transverse cur- 
vature. 

The transverse curvature is discontinuous at E, and 
Eqs. 
quantiti 
terized 
that p 
On the 



12)^( 2.14 ) allow us to compute the shell's surface 
:ies. Because the generators are afhnely parame- 
by —r on both sides of the hypersurface, we have 
= — the shell has a vanishing surface pressure, 
other hand, its surface density is given by 



47]-j'2 ' 



(5.23) 



the ratio of the shell's gravitational mass m to its (ever 
decreasing) surface area Anr'^. Thus far our results are 
virtually identical to those obtained in the first subsec- 
tion. What is new in this context is the presence of a 
surface current j^, whose nonvanishing component is 



3Ma 



J" 



(5.24) 



This comes from the shell's rotation, and the fact that 
the situation is not entirely spherically symmetric. 

To better understand the physical significance of the 
surface current, we express the shell's surface stress- 
energy tensor. 



(5.25) 



in terms of the vector = k" + [j^ / /i) e^. This vector is 
null (in the slow-rotation limit) and has the components 

^"9a " 7 ^* " + 7 (5.26) 



in the coordinates x" = (t, r, 0, 4>) used in A^+; we have 
set 



_ 2Ma 3Ma , 

"fluid = — ^ 1- -X J- 

Zmr 



(5.27) 



The shell's surface stress-energy tensor is now given by 
the simple expression 



(5.28) 



which corresponds to a pressureless fluid of density fi 
moving with a four- velocity £"-. We see that the fluid is 
moving along null curves (not geodesies!) that do not co- 
incide with the shell's null generators. This motion across 
generators is created by a mismatch between fifluid, the 
fluid's angular velocity, and regenerators, the angular veloc- 
ity of the generators. The discrepancy is directly related 
to j-^: 



'relative 



fluid 



■generators 



3Ma 
2mr 



f. (5.29) 



Notice that the fluid rotates faster than the generators, 
which share their angular velocity with zero-angular- 
momentum observers within . Notice also that 
^relative dccrcascs to zero as r approaches 2M: the fluid 
ends up corotating with the generators when the shell 
crosses the black-hole horizon. 

The application considered in this subsection is a varia- 
tion of an analysis presented previously by Musgrave and 
Lake . In some respects our discussion is less general 
than theirs: our black hole is initially nonrotating and it 
does not carry an electric charge. In other respects our 
discussion is more general: here the shell's material is not 
restricted to move along the null generators; as we have 
seen, this generalization leads to interesting phenomena. 



Cosmological phase transition 

In this fourth (and last) application of the formalism, 
we consider an intriguing (but entirely artificial) cosmo- 
logical scenario according to which the universe was ex- 
panding initially in two directions only, but was then 
made to expand isotropically by a sudden explosive event. 
This application is a simplified version of an example pre- 
viously presented by Barrabes and Hogan [ |T9[ . 

The A4~ region of spacetime is the one in which the 
universe is expanding in the x and y directions only. Its 
metric is 



ds^_ = -dt^ + a^{t){dx^ + dy^) + dz 



2 

— 7 



(5.30) 



and the scale factor is assumed to be given by a(t) cx 
t^/^. The cosmological fluid moves with a four- velocity 
= dx"^ /dt, and it has a density and (isotropic) pres- 
sure given by = p- = l/(327rt^), respectively. 

In the region of spacetime, the universe expands 
uniformly in all three directions. Here the metric is 



-dt'^ + a^{t){dx'^ -f dy^ + dz\), 



(5.31) 
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with the same scale factor a{t) as in Ai~ , and the 
cosmological fluid has a density and pressure given by 
p+ = 3p+ = 3/(327rt^), respectively; this corresponds to 
a radiation-dominated universe. 

The history of the explosive event that changes the 
metric from to coincides with a null hypcrsur- 
face E. This surface moves in the positive z± direction, 
and as we shall see, it supports a singular stress-energy 
tensor. The "agent" that alters the course of the uni- 
verse's expansion is therefore a null shell. 

As seen from A4~ , the hypersurface is described by 
t — + constant, and the vector k^da — dt + dz is 
tangent to the null generators, which are parameterized 
by t. Because k'^.^k^ = 0, we have that t is an afRne 
parameter on this side of the hypersurface. The coordi- 
nates X and y are constant on the generators, and we use 
them, together with t, as intrinsic coordinates on E. We 
therefore have — {t,6^), 9^ — {x,y), and the shell's 
induced metric is 



aAB dO^de^ = a^{t){dx'^ + dy^). 



(5.32) 



The remaining basis vectors are e"da = dx, Cyda = dy, 
and Na dx" — —^{dt + dzJ). The nonvanishing compo- 
nents of the transverse curvature are 



C4 B - TT ^AB 



At 



(5.33) 



We note that on the M. side of E, the null generators 
have an expansion given hy 6 = k"^ = 1/t, and that 

Taf3k"kf = P-+P- = 1/{16ttP), where T"'^ is the stress- 
energy tensor of the cosmological fluid. 

As seen from Ai'^, the description of the hypersurface 
is obtained by integrating dt — a(t)dz^, and k'^da — 
dt + a~^dz is tangent to the null generators. We note 
that t is not an afhne parameter on this side of the 
hypersurface: we have that fc^^/c'^ = (2i)-i/c". The 
remaining basis vectors are c^Sq = dx, Cyda = dy, 
Nadx" = —^{dt + adz+), and the nonvanishing com- 
ponents of the transverse curvature are now 



C+-- 
^" " 2t' 



At 



(5.34) 



On this side of E, the generators have an expansion also 
given by 9 = 1/t (since continuity of 9 is implied by 
continuity of the induced metric), and Tapk^k^ = + 
p+ = l/(87rt2). 

The fact that t is an afhne parameter on one side of 
the hypersurface only tells us that the shell must pos- 
sess a surface pressure. In fact, continuity of Cab across 
the hypersurface implies that p is the only nonvanishing 
surface quantity. It is given by 



1 



(5.35) 



the negative sign indicating that this surface quantity 
would be better described as a tension, not a pres- 
sure. The shell's surface stress-energy tensor is S"^ = 
pa^^ e'^e'^. If we select observers comoving with the 
cosmological fluid as our preferred set of observers, then 
—kaU°' = 1 and the full stress-energy tensor of the surface 
layer is given by T^^ = S'"'3(S(i-is), with is denoting the 
time at which a given observer crosses the shell. We see 
that for these observers, —p is the physically-measured 
surface tension. 

Finally, we note that the expressions ~p = l/(167rt), 
9 = l/t, and [Ta^k^k^] = l/ilGnt'^) are compatible 
with the general relation —p9 = [Tapk'^k'^] presented 
in Eq. ([1.5|). This shows that the energy released by the 
shell as it expands is absorbed by the cosmological fluid, 
whose density increases by a factor of /p_ = 3; this 
energy is provided by the shell's surface tension. 
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